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1. Introduction 

We have had very interesting presentations on event-by-event fluctuations during this work- 
shop (see other contributions to these proceedings). This work is based on our recent paper [1] 
and explores from a theoretical viewpoint the correlation data published by the PHENIX [2, 3], 
STAR [4, 5], and CERES [6] Collaborations. These measurements contribute to the previously 
accumulated vast knowledge on event-by-event fluctuations in ultra-relativistic heavy-ion colli- 
sions [6-31]. 

We bring up some basic and rather striking facts manifest in the data of Refs. [2,5] and discuss 
their relevance to the cluster picture of correlations. Let us begin with the PHENIX measurement 
[2] of the event-by-event fluctuations of the transverse momentum at y/s^ =130 GeV. In order to 
simplify the notation, p is used generically to denote the transverse momentum, \pj\ and pi is the 
value of p for the /th particle. Finally, M = Yli^i Pi/n is the average transverse momentum in an 
event of multiplicity n. The statistical treatment of the M variable has been presented in detail in the 
contribution of the speaker to the discussion session (see "Round Table Discussion: Correlations 
and Fluctuations in Nuclear Collisions" in these proceedings [32]). For the reader's convenience 
this material is inluded here in Appendix A. 

2. A look at the data 

The PHENIX results of Ref. [2] are reminded in Table |l|. Several important features of the 
data are immediately seen: firstly, the quantities M, the average in a given centrality class, and 
Op, the inclusive standard deviation of p, are practically constant in the reported "fiducial centrality 
range" c = 0-30%. 

M = const. , Op = const, (at low c). (2. 1) 

This allows us to replace the average momentum at each « by M and the variance at each n by 
Op [I]. For distributions sharply peaked in the n variable we derive in [32] the following formula: 

aS" ^ Op/Vn. (2.2) 

Table 1: Data for the event-by-event fluctuations in the transverse momentum. Upper rows: the PHENIX 
measurement at Jsnn = 130 GeV [2]. Bottom rows: the mixed-event results. 



centrality 


0-5% 


0-10% 


10-20% 


20-30% 


n 


59.6 


53.9 


36.6 


25.0 


On 


10.8 


12.2 


10.2 


7.8 


M 


523 


523 


523 


520 


Op 


290 


290 


290 


289 


Om 


38.6 


41.1 


49.8 


61.1 




523 


523 


523 


520 


—.mix 


37.8 


40.3 


48.8 


60.0 
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Table 2: Verification of the statistical result (2.2) holding for sharp distributions. 



centrality 


0-5% 


0-10% 


10-20% 


20-30% 


—.mix 


37.8 


40.3 


48.8 


60.0 




37.6 


39.5 


47.9 


59.0 



Table 3: Verification of the scaling result (3.2), seen in the approximate constancy of CJdynVn- 



centrality 


0-5% 


0-10% 


10-20% 


20-30% 


O'dynV^ 


60.3 ±1.6 


59.2±1.5 


59.8 ±1.2 


57.7±1.1 



The derivation involves only the statistics and the assumption of sharp distributions, which is suf- 



ficiently well satisfied at RHIC. Corrections to Eq. ( |2.2| ) appearing for broader distributions can be 
obtained. They involve higher moments of n, for instance (y'^{n)/n^. Table ^ shows how well (2.2) 
works - at the level of 1-2%, which is no wonder, just the statistical fact. 

3. The cluster scaling 

Next, we look at the difference of the measured and mixed-event variances of M, which is a 
measure of dynamical fluctuations, 

_2 _ „2 2,mix ,r, , s 

"dyn = '^M~'^M ■ ^) 

Table ^ shows our main observation, namely, that the following scaling holds: 

1 



dyn 



(3.2) 



It works at the accuracy level of 2% in the available centrality range, as can be inferred from 
the approximate constancy of numbers displayed in Table ^ Elementary statistical considerations 
shown in [32] lead to the result 

Pin) " 1 " 

^dyn=L-;^ L COVn{Pi,Pj)^^ X ^^"^ n{p i ^ P j) , (3.3) 

where in the second equality we have used the feature of sharp distributions in n. The quantity 
cov„{pi,pj) is the covariance of momenta pi and pj in events of multiplicity n. 

Now we come to the physics discussion. The scaling ( p^ imposes severe constraints on 
the physics nature of the covariance term. For instance, if all particles were correlated to each 
other, then T!i j=i j=ii^^'^n{PhPj) would be proportional to the number of all pairs, and adyn would 
be independent of n. Thus we see that the combinatorics is truncated - clearly not all pairs are 



correlated, or a finite correlation length develops. A natural explanation of the scaling (3.2) comes 



from the cluster model, depicted in Fig. [T[ The system is assumed to have Nd clusters, each 
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N„| clusters • 




Figure 1: The cluster picture of correlations consistent with scaling (1.4). 



containing (on the average) (r) particles. The particles are correlated if and only if they belong to 
the same cluster, where the average covariance per pair is 2cov*. The number of correlated pairs 
within a cluster is r{r— l)/2. Some particles may be unclustered, hence the ratio of clustered to 
all particles is {Nci){r)/n = a. If all particles are clustered then a = 1. With these assumptions 
Eq. ( p?Tl ) becomes 



■'dyn 



a{r{r- 1)) 
{r)n 



ar 



-cov 



-cov 



(3.4) 



where we have introduced r* = {r{r— l))/(r), the ratio of the average number of pairs in the cluster 
to the average multiplicity of the cluster. For a fixed number of particles in each cluster we have 
r* = (r) — 1, for the Poisson distribution r* = (r), while for wider distributions r* > (r). Equation 



( |3.4| ) complies to the scaling ( |3.2| ) as long as the product ar*cov* does not depend on n (in the 
fiducial centrality range). This is the basic physics constraint that follows from the data. For the 



reasons discussed above it is appropriate to term the scaling (3.2) the cluster scaling 



4. Cluster scaling with other correlation measures 



There are numerous quantities used as measure of the event-by-event fluctuations. Fortunately, 
as shown in [32], for sharp distributions and for the case where the dynamical correlations are small 
compared to the statistical correlations, adyn ^ Om, all the popular measures are simply related to 
each other, or to the covariance, which is the basic quantity designed to describe correlations. We 
have 



Pt 



F = 



^dyn 



_2 Ly;COV7r(/7;,/7y) 
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p^n^ 
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L-^;COVn(p;,P;) 

InOp ' 

(ApApj) = Y,coyn{Pi,Pj), (4.1) 

where in the definition of O^j, one uses Sn = pi + . . . + p,,. The last equality states the fact that 
the STAR measure ApiApj is just the statistical estimator of the covariance, see [1]. The cluster 
scaling thus manifests itself for the various measures in the following manner: 



y 






\fn 


^Pt 


~ 1, 


^PT 


~ 1, 


{ApiApj) 


~ 1. 



(4.2) 

Since this is an important point, we stress again: since at RHIC the distributions are 1) sufficiently 
sharply peaked in n and 2) the dynamical fluctuations are small compared to the statistical fluctua- 
tions, all popular measures of event-by-event fluctuations are all proportional to each other and to 
the covariance. One may pass from one measure to another without difficulty. If conditions 1) and 
2) are not satisfied, as may be the case at lower energies, periferal collisions, etc., then of course the 
measures remain different and some may be better tuned for certain analyses. Full information on 
correlations could be acquired by simply evaluating the covariance Y^i^j cov„ separately for each 
n. Even if 1) or 2) are relaxed, the correlation measures are still related to the sum of the weighted 
covariances at various n, with the weights dependent on the particular measure. The necessary 
formulas can be very straightforwardly derived along the lines of [32]. We urge that such a study 
of the dependence of covariance on n be made on data with sufficiently large samples. 



5. Cluster scaling in various experiments 

Now let us have a look at various experiments. The PHENIX results at y^S^w =130 GeV [2] 
discussed in Sect. || comply nicely to the cluster scaling. The same is true of the STAR data (see 
Fig. 3 of [5]) at y^Siwv = 200 GeV and 62 GeV, where the quantity dn/dr] (ApiApj) flattens out for 
large centralities. For ^/snn = 130 GeV and 20 GeV the STAR data seems to somewhat depart from 
the scaling, with sizeable error bars for 130 GeV. The PHENIX data at ^/s^n = 200 GeV shows 
non-monotonic behavior [3]. The STAR data at ^/snn = 200 GeV also shows non-monotonic 
behavior at large centralities, but only for the analysis where all produced particles are taken into 
account. When the analysis is constrained to the in-plane or out-of-plane particles, then the cluster 
scaling is satisfied remarkably well (see Paul Sorensen's talk showing the STAR preliminary data). 
This is a very intriguing phenomenon that needs to be understood. The CERES data [15], using 
the Epj- variable, also complies to the cluster scaling within the error bars. Recall that according to 
( p!^ ) we look for Z^j, ~ 1/ In conclusion, the cluster scaling is seen in some measurements to a 
remarkable accuracy, while is other non-monotonous behaviour is apparent. The situation requires 
careful clarification. 
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6. Strength of correlations 

Before performing the analysis of clusters in a more quantitative manner we need to consider 
the effects of acceptance and the detector efficiency. This is particularly important in the event- 
by-event analysis, since the experiments select particles with very clearly identified tracks, hence 
the detector efficiency a is reduced. The number of observed particles is proportional to a, and 



the number of pairs contributing to the covariance is proportional to a^. Thus Eq. (3^) may be 
rewritten as 

r* r* 

^lyn = — COV*=a— COV*, (6.1) 

' n full n obs 

where "full" denotes all particles (that would be observed with 100% efficiency), while "obs" 
stands for the actually observed multiplicity of particles. Thus 

'^OV =^dyn^- (6-2) 

Our estimate for a in the PHENIX experiment is of the order of 10%, which together with the 
numbers of Table [l] gives 

, 0.035 GeV^ 
COV ~ . (6.3) 

In the considered problem the coefficient 0.035 GeV^ is not a small number when compared to the 
natural scale set by the variance Op ~ 0.08 GeV^ (we recall that | cov* | < a^). 

Very similar quantitative conclusions are reached with the STAR data. Taking the values from 
Table I of Ref. [5] and guessing a = 0.75 we find cov*r* = 0.058,0.043,0.035,0.014 GeV^ for 
^/s^ = 200,130,62 and 20 GeV, respectively. The value at 130 GeV is close to the PHENIX 



value (5.3). Interestingly, we note a significant beam-energy dependence, with cov*r* increasing 
with ^snn. This may be due to the increase of the covariance per correlated pair with the increasing 
energy, and/or an increase of the number of particles within a cluster. 

Using the above numbers we find that for r* = 1 (for instance the case where all clusters have 
just two particles) the value of cov* assumes almost a half of the maximum possible value, a^. This 
is very unlikely, as dynamical estimates presented below give cov* of the order at most 0.01 GeV^. 



Thus a natural explanation of the values in ( |6.3| ) is to take a significantly larger value of r* - just 
put more particles inside the cluster. Of course, the higher value, the easier it is to satisfy ( |6.3[ ) even 
with small values of cov* . This scenario will be elaborated in the next section. 



7. The nature of clusters 

Several mechanisms can be brought up to describe the formation of clusters in the momentum 
space. The basic physics question concerns the nature of clusters. Here we discuss a few popular 
scenarios. 

7.1 Jets 

The most prominently explored mechanism is the formation of (mini)jets. In Ref. [1] we have 
discussed this issue, with the conclusion that the explanation of the centrality dependence of the pj 
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Figure 2: The covariance cov*gj. of the two pions produced in the decay of the resonance om mass M. 



fluctuations in terms of jets based solely on scaling arguments is not conclusive. Any mechanism 
leading to clusters in pT would do. Microscopic realistic estimates of the average number of jets 
and the magnitude of the covariance of particles originating from a jet are necessary here, including 
the interplay of jets and the medium. For the current status of this program the user is referred 
to [31,33]. 

7.2 Resonance decays 

Another natural mechanism for momentum correlations is provided by resonance decays. Re- 
call that resonances are crucial behind the success of the statistical approach to particle production 
in relativistic heavy-ion collisions. We have made an estimate of this effect in the thermal model 
of Ref. [34, 35]. The covariance cov*gj, of the two pions produced in a resonance decay is given by 
the formula 



jd^pj^IPij^lPi5(^){p-p,-P2)c'-^ 



(7.1) 



where dNR/d^p is the resonance distribution in the momentum space (obtained from the Cooper- 
Frye formula as described in Ref. [36]), M is the mass of the resonance, pi and p2 are the momenta 
of the emitted particles, Ep is the energy of a particle with momentum p, and the function C 
represents the appropriate experimental cuts. Here we take C as for the PHENIX experiment, 
with 0.2 GeV < /jfj < 1-5 GeV, and |3;| < 0.35. The average momentum of the pion is {p"^) = 
523 MeV. The parameters of the model are T = 165 MeV (the universal freeze-out temperature), 
Pmax = 7.15 fm (the transverse size of the system), and t = 7.86 fm (the proper time at freeze- 
out). The numerical results presented in Fig. |^ show that for the resonance mass between 500 MeV 
and 1.2 GeV the covariance cov*^^ varies between 0.01 GeV^ at low masses to —0.018 GeV^ 
at M ~ 1.1 GeV. These numbers are small compared to the scale ~ 0.08 GeV^. Moreover, 
there is a change of sign around M = 700 MeV and cancellations between contributions of various 
resonances are possible. In fact, a full-fledged simulation with Therminator [37] revealed a 
completely negligible contribution of resonances to the pj correlations with model parameters 
appropriate for the STAR and PHENIX experiments. 
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Table 4: Predictions of the "lumped cluster" model. The quantity T is the universal freeze-out temperature, 
(jSr) is the average transverse flow velocity, cov*ji is the average covariance of the pion pair, and r* is 
an estimate for the average number of charged particles in the cluster based on the PHENIX results for 
130 GeV. 



T [MeV] 


10 


100 


120 


140 


165 


200 


(i8r> 


0.94 


0.72 


0.69 


0.58 


0.49 


0.31 


oj [GeV2] 


0.056 


0.19 


0.15 


0.15 


0.14 


0.12 


cov*^ [GeV2] 


0.027 


0.011 


0.0088 


0.0063 


0.0034 


0.0006 


J,* _ 0.035 GeV^ 


1.3 


3.2 


4.0 


5.6 


10.3 


58.3 



7.3 Lumps of matter 

Finally, let us consider a model of momentum correlations which assumes that the particle 
emission at the lowest scales occurs from local thermalized sources. We call this picture the 
"lumped clusters": lumps of matter move at some collective velocities, correlating the momenta of 
particles belonging to the same cluster, see Fig. |l|. Each element of the fluid moves with its col- 
lective velocity and emits particles with locally thermalized spectra. This picture was put forward 
as a mechanism creating correlations in the charge balance function [36, 38] resulting from charge 
conservation within a local source. The covariance between particles / and j emitted from a cluster 
moving with a velocity u is given by the equation 

IdL,U^Jd'p,{pJ-{p^))fl'{pi)Id'p2{pl-{p'))fnP2) 
IdZ,U^!d^p,fnPl)Id'P2f}'iP2) ' 

where f"{p) = (exp(;? • u/T) it 1)^' is the thermal distribution in the local reference frame and 
dL^ denotes integration over the freeze-out hypersurface. In this calculation we adjust the average 
transverse flow velocity (jSr) at each T such that the slope of the pion spectra agrees with the data. 
Of course, lower t requires higher (jSj-). The result turns out to depend strongly on the temperature. 
For the emission of correlated pion pairs one gets the results shown in Table ^ The last row of 



the table gives the number of particles based on formula ( |6.3| ), obtained for PHENIX at ^/sf^ = 
130 GeV. We note that for realistic values of thermal freeze-out parameters the experimentally 
estimated value of the covariance cannot be accounted for, unless the number of charged particles 
belonging to the same cluster is at least of the order 4—10 (assuming the Poisson distribution). For 
wider distributions in the variable r a lower number is requested. Thus the number of all particles 
(charged and neutral) belonging to a cluster is estimated as 

4-6-15, (7.3) 

which is one of our main results. 

8. Dependence on maximum transverse momentum 

An interesting result is obtained when the upper limit of integration in the transverse momen- 
tum, p™^", is imposed. Figure || shows the dependence of Fpj. on pf'^'^ and compares the result to 
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• •• 
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2 




1 



•0.5 



1 



1.5 



2 



2.5 



3 



Pt 



[GeV] 



Figure 3: Dependence of Fp^ on the cut in the transverse momentum. The data points come from Ref. [3]. 
Red points show the model calculation with T = 165 MeV, while the black points are for T = 130 MeV. 

the PHENIX data. The red points show the model calculation with T = 165 MeV and ar* =2.1, 
while the black points are for T = 130 MeV and ar* = 2.9. The value of ar* is tuned such that the 
model describes the data well. For realistic values of the acceptance a we get values of r* of the 
same order as in estimate ([7J|). 

A feature of our model is a saturation of Fpj, at large values of pj^'^. This saturation is a con- 
sequence of the thermal model. As a matter of fact, all relevant quantities saturate when integrated 
over pj. Unfortunately, the experimental data do not extend to the saturation region. What is some- 
what surprising at the first sight is the rather large value of p™^" needed for saturation: about 2 GeV 
for r = 130 MeV and about 3 GeV for T = 165 MeV. Figure ^ shows the results for the average 
momentum, its variance, and the covariance resulting from the lumped cluster model. We note that 
all these quantities saturate at rather large values of /j™^". In that sense in the present context the 
"soft" thermal physics reaches transverse momenta up to 2-3 GeV! 

9. Conclusion 

We list our main results: 

1. The cluster scaling of ajyj, with l/nin the fiducial centrality range 0-30% at PHENIX points 
at the cluster picture of the fireball. 

2. Similar scaling can be also seen at STAR [4,5] (see also Paul Sorensen's talk) and at CERES 



3. ScaUng can be found with various correlation measures, which in the limit of sharp distribu- 
tions and small dynamical compared to statistical fluctuations become equivalent and related 
in simple ways to the covariance, as discussed in [32]. 

4. The clusters may a priori originate from very different physics: (mini)jets, droplets of fluid 
formed in the explosive scenario of the collision, or other mechanisms leading to multiparti- 
cle correlations. 



[6]. 
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0.12 
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0.004r 
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0.001 

^[GeV] 
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Figure 4: The dependence of the average momentum (top), its variance (middle), and the covariance (bot- 
tom) resuhing from the lumped cluster model for T = 130 MeV (left) and T = 165 MeV (right) on the 
cut-off in the transverse momentum, p™"". 



5. (Mini)jets just produce clusters, so it is impossible to prove or disprove their existence based 
solely on the centrality dependence of the correlation data at soft/medium values of pr- 

6. Resonance decays in a thermal model yield a very small value of the pr-covariance. 

7. In a model where matter forms lumps moving at similar collective velocity and particles 
move thermally within a cluster, description of the RHIC data requires on the average 6-15 
particles in the cluster In general, a larger number of particles within a cluster helps to obtain 
the large (compared to Op)/h) measured value of o^^^. 

8. In the "lumped cluster" model the p^'^^ dependence of Fp-^- grows monotonically and then 
saturates. It would be interesting to confront this finding to the data. 

9. Detailed microscopic modeling would be very useful in order to better understand the prob- 
lem on transverse momentum correlations in relativistic heavy-ion colUsions. 
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A. Relation of measures of fluctuations to covariance 

One difficulty in comparing results of event-by-event momentum fluctuations presented by 
various experimental groups is the multitude of measures used. Here we briefly show how under 
assumptions of 1) small dynamical compared to statistical fluctuations and 2) sharp distribution in 
the multiphcity variable these measures are simply proportional to the covariance. Although these 
remarks are perhaps obvious to practitioners in the field, they seem worth reminding, as discussions 
showed confusion. 

Suppose we have events of class n (formally, this can be any number chracteristic of the event: 
the multiplicity of detected particles, the number of participants, the response of a given detector, 
etc., distributed according to the probability distribution Let p„{pi,..., p„) denote the n-particle 
distribution of variables pi within events of class n, e.g, the distribution of transverse momenta in 
events of a fixed multiplicity n). The subscript „ indicates that p depends functionally on n. The 
full probabihty distribution of obtaining event of class n with momenta , . . . ,p„ is 

PnPn{Pl,---,Pn)- (A.l) 

The marginal probability distributions are obtained from Pn{p\,... ,Pn) by integrating over k mo- 
menta, 

PniPU- ■ ■ ,Pn-k) = J dpn-k+\ ■ ■■dpnPni.PU- ■ ■ ,Pn)- (A.2) 

Next, we introduce the relevant moments for the distributions of class n: 

P^ = Jdpp„{p)p, 

<^niP) = J dppn{p){p -P^f , 

coy„{pi,p2) = J dpidp2Pnipi,P2)ipi-p^)ip2-p^), (A.3) 

where Pn{p) and pniPijPi) are the one- and two-particle marginal distributions within the class n. 
Now, in a typical setup we are interested in broader classes, containing n in the range n\<n< n2. 
We denote for brevity = . 

Let us illustrate the basic statistical facts on the example of the measures (7jy„ and Fpj^ . For 
other measures the analysis is analogous. Consider the variable M„ = (pi Pn)/n, i.e. the 

average value of the variable p. Then 

M = Y^Pn I dpi... dpnPn{p\ ,---,Pn) = Y^^nP"^ 

r \ n 

= Y,Pn / dpi...dpnPniPl,---,Pn)— Y [(Pi -P^)iPj -P^) +Pn^] ' 
n J " ,j=l 

al,=W-M^=YPn-p;.^-(YPnT^ +YPn^^ + "LPn^ L COVn{pi,p j) . {A.4) 
n \ n J n " iW=l 

Suppose mixing of events is performed. Then, by definition, no correlations are present, i.e., 
coN'^{pupj)=0, and 

0lr'^ = YPnP^'-{^PnP-^^ + LPn^. (A.5) 
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By definition, o^^^ = o^ — CJ^™"- Witli above results 

f^dyn = L ^« ~2 L '^°^« 'Pj)- ^) 

Note that all above results are exact, just following from obvious manipulations. Thus o^^^ is 
a weighted sum of total (summed over particle pairs) covariances at fixed n, Y.i^jCOWniPi,Pj)> 
with weights equal to Pn-^- Moreover, all quantities in Eq.(A.4) or i ^.5\ > are possible to obtain 
experimentally from the given event sample. 

Now let us have a look on Fp^^ = (a/w — -y/oVx ) / -y/oVx , where (O = o^/M. We find imme- 
diately 



1 + (A.7) 

2,mix ^ ' 



Again, this is an exact relation. At RHIC a|y^ <^ '^m™"' hence we can expand 

Fp^^l-^, (A.8) 

which shows the proportionality of the two measures in the limit of small dynamical correlations. 

A further simplification occurs when the distributions P„ are sharply peaked around some n, 
which again is sufficiently well satisfied at RHIC. Then for a smooth function f{n) 



(A.9) 



n " 



In this sharp limit 



^2M. ~ (A. 10) 
n 

For other measures the results are similar. For we have from definition S^^, = o\^J'p'^, for 
Opj, under conditions 1) and 2) 

^P, ^ \/ ^ - <y^^P) ^ 'I' I (A.11) 
V " ^nOn(p) 

where Sn=p\ + ..- + Pn- 

Conclusion: Since at RHIC conditions 1) and 2) liold, tlie popular measures of event-by- 
event fluctuations are proportional to tlie covariance. Full information on correlations could 
be acquired by simply evaluating the covariance L,yjCov„ for each n. If 1) or 2) are relaxed, 
then, of course, the measures are no longer equivalent, but they are still related to the sum of the 
weighted covariances at various n in the way dependent on the particular measure. 
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Finally, we make a digression concerning the inclusive distributions, not used in our deriva- 
tions but appearing frequently in similar studies. They should not be confused with the marginal 
distributions, to which they are related as follows: 

Pin(p) = Y^Pn / dpi . . .dpnY,S{P - Pi)Pn{PU ■ ■ ■ ,Pn), 
n -I !=1 

/n 
dpi...dpn 5{p-pi)5(q-pj)pn{pi,...,pn), (A.12) 

with the properties 

J dppinip) = n, 
J dpdqpm{p,q) = n(n-l), 

p- = Y,nPnP^- (A. 13) 

n 
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